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ABSTRACT 


We present our results of a combined analysis of radial velocity and light 
curves of a double-lined spectroscopic and eclipsing binary KIC 09246715, ob¬ 
served photometrically by the Kepler satellite, and spectroscopically with the 
OAO-1.88m telescope with the HIgh-Dispertion Echelle Spectrograph (HIDES). 
The target was claimed to be composed of two red giants, one of which is showing 
solar-like oscillations. We have found that the mass and radius of the primary 
are Mi = 2.169 ± 0.024 Mq and Ri = 8.47 ± 0.13 Rq, and of the secondary: 
M 2 = 2.143 ± 0.025 Mq and R 2 = 8.18 ± 0.09 Rq, which conhrms the double¬ 
giant status. Our secondary is the star to which the oscillations were attributed. 
Results of its previous asteroseismic analysis are in agreement with ours, only 
signihcantly less precise, but the subsequent light-curve-based study failed to de¬ 
rive correct mass and radius of our primary. KIC 09246715 is one of the rare 
cases where asteroseismic parameters of a solar-like oscillator were conhrmed by 
an independent method, and only the third example of a Galactic double-giant 
eclipsing binary with masses and radii measured with precision below 2%. 

Subject headings: binaries: eclipsing — binaries: spectroscopic — stars: evolution — 
stars: fundamental parameters — stars: oscillations (including pulsations) 
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Introduction 


The unprecedented photometric precision of the Kepler mission flBorucki et alJl2010l) 
has opened new possibilities in various helds of astronomy, including studies of eclipsing 
binaries and asteroseismology. The former has been known for many decades to provide 
accurate and precise stellar parameters, sometimes not possible to derive with other 


meth ods (e.g. absolute radii). Since the discovery of pulsatin g sdB stars flKilkennv et ah 


1994t) and solar-like oscil 


ations in stars other than the Sun (IBrown fc Gilliland 


1994 


Kieldsen fc Bebbindll995fl . asteroseismology became another technique capable of providing 


such results. Until the launch of CoRoT and Kepler satellites the number of other stars 
showing the solar-like oscillations was very low, but now it is counted in thousands, and 
contains main sequence, as well as giant stars. 

Absolute values of fundamental stellar parameters have multiple applications 
i n modern astronomy, but they need to be known with precision of at least 2-3% 


flLastennet fc Vails-Gabaud 


2 OO 2 I) . The best possible characterisation comes from double- 


lined spectroscopic and eclipsing binaries, that also show pulsations. Such systems are 
extremely useful for testing stellar evolution models, as the masses and radii can be, 
in principle, measured with two independent techn iques. Around 50 ec lipsing systems 


containing oscillating giants have been reported by 


Gaulmeetah 


( 20131. hereafter: G13), 


but without spectroscopic follow up. Only KIC 08410637, c omposed of an oscillat ing giant 
and a main-sequence star, has been studied in more detail ( Frandsen et ah 20131 . and a 


very good precision in absolute parameters has been achieved. Despite being very useful in 
general, currently the asteroseismic measurements usually do not give results of sufficient 
precision. Nevertheless, the large number of known oscillating stars undoubtedly shows the 
importance of this technique. Discovery and characterisation of the “keystone” double-lined 
spectroscopic, eclipsing and oscillating systems will help to test and enhance the capabilities 
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of asteroseismology directly. In this work we present such an example. 


2. The target 


The system KIC 09246715 (HD 190585, V = 9.65, hereafter: K0924) was identihed 
as an eclipsing binary by the Kepler miss i on, and was listed i n the Kepler Eclipsing 


Binaries Catalog (KEBC; 


Prsa et ah 


2011 


Slawson et ah 


201lh . It was not reported as 


a double-lined spectroscopic binary till now. According to G13 and 


Gaulme et al 


1I2OI4J. 


hereafter: G14), it is composed of two red giants, one of which is showing solar-like 
oscillations with peak frequency and separation of I'rnax = 102.2 /xHz, Au = 8.3 /xHz (G13), 
or Umax = 106.38(75) /iHz, Au = 8.327(10) /xHz (G14). Activity is also reported, in a form 
of periodic {Pyar = 93.3 d) brightness modulations coming from spots. 


Using asteroseismology G13 and G14 Erst derived mass and radius of the oscillating 
star. Then they combined it with results of light curve modelling (i.e. Ri/a, R 2 /a) and 
Kepler’s 3rd law, inferring the properties of the other component: 2.06(13) Mq, 8.10(18) Rq 
for the oscillating component (their primary), and 1.1(3) Mq, 6.8(1) Rq or 3.3(5) Mq, 
9.7(2) R 0 for the other, but described the latter solution as less probable. With our 
spectroscopy, we are able to directly revise their results, confronting the asteroseismology 
with direct modelling of the light and radial velocity (RV) curves of this system. 


Eclipsing binaries composed of two red giants, like K0924, are relatively rare, and 


not many such systems have t 


The on-line DEBGat catalog (jSouthworth 


leir components 


parameters measured with high precision. 


20151 ). which lists detached eclipsing binaries 


with masses and radii measured do wn to 2-3%, contains only t wo such Galactic systein s 


HD 187669 f|Helminiak et al. 


2015 ) and ASAS J1800 52-2333.8 (Suchomska e 


and a handful of others found in LMG and SMG (e.g. 


Pietrzvhski et al. 


2013 


: al 


20151) - 


Graczvk et al. 
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2014f) . Therefore K0924 is also important for studying the late stages of stellar life and 


testing the models of evolution and structure. 


3. Data and methodology 

The system K0924 is one of the targets we have observed as a part of a larger 
spectroscopic project, aimed for characterisation of northern detached eclipsing binaries, 
including objects selected from the Kepler held. All of them have been observed and 
analysed in a similar way, and the detailed description of the observing procedure, data 
reduction, stability tests, orbit and light curve htting, etc. will be explained in a forthcoming 
paper (Helminiak et ah, in prep.). However, we shortly present them below. 

We follow the convention that the primary component is the one being eclipsed during 
the deeper minimum (the primary eclipse). In our case, the primary is the slightly more 
massive and larger star. 


3.1. Observational data 


In this study we make use of the publicly available Kepler mission photometry. We 
used the de-trended relative hux measurements fdtr, that were later transformed into 
magnitude difference Am = —2.5\og{fdtr), and hnally the KEBC value of kmag was added. 
The de-trended data were downloaded directly from KEBC, where currently only the 
long-cadence measurements for this star are available. 


A total of 8 spectra were taken during several runs between July 2014 and late 
April 2015, at the 1.88-m telescope of the Okay ama Astrophys ical Observatory with the 


High-Dispersion Echelle Spectrograph (HIDES; 


Izumiura 


19991) . The instrument was fed 
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through a circular hber, for which the light is collected via a circular apertur e of projected 


on-sk y diameter of 2.7 arcsec, drilled in a flat mirror that is used for guiding flKambe et ah 


20131) . An image sheer is used in order to reach both high resolution {R ~ 50000) and good 


efficiency of the system. Wavelength calibration was based on ThAr lamp exposures taken 
every 1-2 hours, which gave the stability of the instrument at the leyel of 40-50 m s“^, as 
measured from multiple observations of four RV standards. Each science exposure was 1500 
seconds long, but due to unstable atmospheric conditions, the resulting signal-to-noise ratio 
(SNR) usually varied between 75 and 105 (63 at one time). Spectra were reduced with 
dedicated IRAF-based scripts, that cope with the mosaic character of the HIDES detector. 
Reduction included correction for overscan, bias, flat field, cosmic rays and bad pixels, as 
well as extraction and wayelength calibration of 53 orders (4360-7535 A). 


3.2. Light curve analysis 

For the light curve (LC) fit we used version 28 (v28) of the code JKTEBOP 


fSouthworth et ah 

2004a, 

3 ), which is based on the EBOP program ( 

Popper & Etzel 

1981) 


The complete long-cadence Q0-Q17 curve was first used to find the best-fitting parameters, 
including period P, primary (deeper) eclipse mid-time Tq, eccentricity e, periastron 
longitude u, inclination i, ratio of fluxes L 2 /L 1 , ratio of central surface brightnesses J, sum 
of the fractional radii ri -|- r 2 (in units of major semi-axis a), and their ratio k. A small but 
clear increase in brightness was noticed around both eclipses, so we also fitted reflection 


coefficient s, and got 7.4(1.0) x 10 ^ anc 


6. 3(1.0) X 10 Logarithmic limb darkening 


(LD) law (IKlinglesmith fc Sobieskilll970l ) was assumed, with coefficients (fixed in the fit) 


interpolated from the the tables published on the PHOEBE websitqj- For this, the gravities 


^http://phoebe-project.org/1.0/?q=node/110 
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log(5') were taken from an initial fit (corrected la ter), which results, i.e. masses and radii, 


were then compared to the PARSEC isochrones flBressan et ah 


2 OI 2 I) in order to estimate 


the temperatures. The gravity darkening coefficients were always kept hxed at the values 
appropriate for stars with convective envelopes {g = 0.32). At the end, the code calculates 
the fractional radii ri_ 2 , and fluxes Li, 2 - As input we also used spectroscopic flux ratios as 
found by TODCOR (see Sec. 13.31) . 

The JKTEBOP does not £t for spots or oscillations, however, it offers a number of 
algorithms to properly include systematics in the error budget. We assumed that after 
years of nearly continuous observations they are averaged out over the orbital period, and 
decided to treat them as a corr elated noise and run th e residual-shifts (RS) procedure to 


calculate reliable uncertainties (jSouthworth et ah 


2011 ). To run RS on the whole Q0-Q17 


curve would take almost two weeks on the computer we used, so for the errors estimation 
we decided to split the data and analyse each chunk separately, and then calculate weighted 
averages of the resulting parameters. Due to its long period, reaching almost the time 
span of two quarters of Kepler data, both eclipses were not always visible during the same 
quarter. Adjacent quarters with only primary and only secondary minimum were merged 
together into one set and cropped, so the resulting light curve time span was around 90 days 
and covered both minima. This was done for Q6-I-Q7, Q8-I-Q9, QlO-l-Qll, and Q12-I-Q13. 
Both eclipses were observed in Q3, Q5 and Q14. No eclipses were seen in QO, Ql, Q2, Q4, 
Q15, Q16 and Q17, which were rejected from the RS stage. Because of the long period, 
P,To,e,oj and reflection coefficients, as well as the LD coefficients, were hrst held hxed, 
only perturbed later during the proper RS stage. The hnal parameter errors were obtained 
by adding in quadrature the formal error of the weighted average and the rms of the results 
for each quarter. 









3.3. Radial velocities and orbital fit 


RVs were measured with o ur own implementation of the two-spectra cross-correlation 


function technique TODCOR fIZucker fc Mazeh 


with ATLAS9 and ATLAS12 codes ( Kurucz 


R)94) with synthetic spectra computed 


I 992 I) as templates. Single measurement 


errors were calculated with a bootstrap approach (IHelminiak et ah 


[a m 


Zucker fc Mazeh 


2012 ). Flux ratios 


199J) were also calculated, and used as input for JKTEBOP. For 


this, only orders from the Kepler bandpass were used, with weights corresponding to the 
total response factor at a wavelength equal the center of a given spectral order. The RVs 
obtained from our HIDES spectra, together with their errors and SNR of the spectra, are 
listed in Table [H 


EDITOR: PLACE TABLE [U HERE. 

The RV measurements were analysed with the code V2FIT, which hts a double- 
Keplerian orbit by using the Levenberg-Marquartd algorithm. As free parameters, we 
set the time of periastron passage Tp, the velocity semi-amplitudes K 12 , the primary’s 
systemic velocity 71 , and the difference in components systemic velocities 72 — 71 . Period 
P, eccentricity e and periastron longitude u were kept hxed on values found in JKTEBOP. 
In order to hnd reliable uncertainties, that include the influence of systematics coming 
possibly from the pulsations and low number of spectra, we run 10000 bootstrap iterations. 
As expected, we found the systematics to be the dominant source of parameter errors. 


3.4. Absolute values of parameters 

The absolute values of parameters and their uncertainties were calculated with 
the JKTABSDIM code, available together with JKTEBOP. This simple code combines 
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the spectroscopic and light curve solutions to derive a set of stellar absolute dimensions 
(Ml 2 , Ri, 2 i a), and related quantities (Li 2 /-^©, log( 5 'i, 2 )). If desired, it also calculates 

distance, but requires multi-color photometry, E{B — V), and both effective temperatures 
as input. Due to lack of such data, we did not attempt to estimate the distance. 


4. Results 

The light and RV curves are presented in Figures [T] and [H respectively. Values of 
stellar parameters can be found in Table [2l We have found that the mass and radius of 
the primary are Mi = 2.169 ± 0.024 M© and Ri = 8.47 ± 0.13 R©, and of the secondary: 
M 2 = 2.143 ± 0.025 M© and R 2 = 8.18 ± 0.09 R©. K0924 is composed of two very similar 
stars that evolved significantly from the main sequence. We reached a good precision 
of ~1.1% in masses and 1.1-1.5% in radii, which makes our results valuable for further 
comparison with models of stellar structure and evolution, and places K0924 in a very 
small group of double-giant systems with precisely measured masses and radii. Precision in 
masses is mainly hampered by the low number of spectra, and likely by the influence of 
spots, oscillations, and instrument distortions on the RV measurements, which may explain 
the apparently non-random residuals. Precision in radii is lowered by the oscillations, which 
we did not model and reduce, but treated as a red noise in the light curve, and included in 
the error budget. 


EDITOR: PLACE TABLE |2] HERE. 


EDITOR: PLACE FIGURE □ HERE. 


EDITOR: PLACE FIGURE [2] HERE. 
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These results are different from G13, who give 1.7(3) Mq and 7.7(4) Rq for the 
oscillating component, and 0.8(7) Mq and 5.9(3) Rq for the other, as well as from G14, 
who give, analogously, 2.06(13) Mq + 8.10(18) Rq, and 1.1(3) Mq + 6.8(1) Rq. Note that 
both previous studies adopt the opposite definition of primary/secondary, and solar-like 
oscillations are found on the star that we define as the secondary, i.e. slightly less massive 
according to our orbital solution. First, both stars have almost equal masses and radii, 
while the solution of G14 gives M 2 /Mi ~ 1.87 and k ~ 1.19. Our mass ratio close to 1 
comes directly from the spectroscopy (Fig. [2]), and its value is very robust. The ratio of 
radii was constrained using the spectroscopic flux ratios from TODGOR, which we also find 
reliable - a quick inspection of any of the spectra reveals that the lines of both components 
are of similar depths, meaning similar effective temperatures and levels of continuum, and 
the peaks of cross-correlation function were of basically the same height. Finally, a quick 
comparison with evolutionary models (see next Section) shows that there is no isochrone 
that reproduces the results of previous studies. This means that the asteroseismology 
provides correct, yet less precise, stellar parameters, but light curve solutions should be 
supported by spectroscop ic values of ma s s and flux ratios. It has also been found in the 
case of KIG 08410637 by 


Frandsen et ah 


fl2013l) . Moreover, some conclusions presented in 


G14 should now be treated with caution. 


It is worth to note that despite very large separation, the system might have reached 
or is close to a form of a tidal equilibrium. If the reported period P^ar = 93.3 d is related 
to the rotation, it would coincide with the value of Prot,ps = 94.2(3) d, which is the 


value of rotation expect ed for the 


pseudo-synchronisation (iHut 


1981 


given orbita 


Mazeh 


2008 


period and eccentricity in the state of 
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Discussion 


In Figure [3] we compare our results , i.e. masses and radii, with the theoretical 


PARSEC isochrones flBressan et al 


20121) . that include calculation of absolute magnitudes 


in the Kepler photometric band. In this set of models the solar metallicity is reached for 
Z=0.0152. 


EDITOR: PLACE FIGURE [3] HERE. 


From Figure |3] one can deduce that the system may be ~950 Myr old, if solar 


metallicity is assumed. At this stage, however, t 


on the mass-radius plane flHelminiak et al 


le age and [M/H] are strongly degenerated 


20151) . Both components are probably on the 


red giant branch, before the red clump stage. It would be possible to reproduce both the 
observed radii on the red clump, by assuming a lower metallicity {[M/H] < —1.6), but 
then the system would be even younger (600-700 Myr) and of much earlier spectral type 
than suggested by the KEBC effective temperature of 4699 K. Both components likely 
have nearly equal Te//’s, and the spectra show many features, so we find the ~4700 K, 
~950 Myr and solar metallicity solution more probable, however with some caution about 
the temperature itself (see further text). 

We also checked the value of metallicity of —0.39 dex, given in the Mikulski Archive for 
Space Telescopes (MAST)I^, which was, however, obtained from analysis of combined light 
of both stars. It is listed together with log( 5 ') = 2.42, which is not in agreement with our 
results, so we treat this value of [M/H] with caution, but we cannot exclude it completely. 
The best-fitting PARSEC isochrone for this [M/H] is for the age of 780 Myr, and is also 


^http://archive.stsci.edu/index.html 
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shown in Fig. |3l For completeness, we also present an isochrone for a more metal abnndant 
case [M/Ff]=0.20. The best-fitting age is then 1.04 Gyr 


Dne to rapid changes in the stellar strnctnre dnring the red giant phase (compared to 
the main seqnence), the isochrones predict a variety of possible temperatnres within the 
range of masses we fonnd for K0924. However, the evolntion in T^ff is related to the growth 
of the star, so we can nse the radii estimates to predict temperatures of both components. 
However, due to lack of any [M/H] estimate that we find reliable, we can not pick any 
certain scenario. On the 950 Myr, solar-metallicity isochrone, we find = 4990(11) K, 

and = 5013(7) K. On the 780 Myr, —0.39 dex isochrone we get = 5250(20) K, 

and Tf,ff ^2 = 5300(12) K. Finally, the 1.04 Gyr, 0.20 dex gives = 4876(9) K, and 
Teff ,2 = 4896(7) K. To obtain temperatures around 4700 K, one would have to use an 
isochrone for even higher metallicity. The difference in predicted temperatures is large 
enough for the modern spectral analysis methods to distinguish between solar and sub-solar 
cases. Note also that the ratio of predicted yalues of temperatures is much closer to 1 


than 1.04, which was found by (iGaulme et ah 


2014J ) in the light-curve fit, and that in their 


solution the oscillator (our secondary) is cooler. The K0924 system would benefit from 
further spectroscopy (to disentangle the component spectra) and multi-band photometry, 
that altogether would allow for precise temperature and metallicity determination. It would 
also improve the age estimation, and precision in masses and radii even more. 


Finally, taking into account the similarities of the two components, we suspect that 
both stars may show similar oscillations that blend in the light curve, and may have 
confused G14, although the mass and radius of their pulsating component, and the 2014 
values of e,i and (i?i -|- i? 2 )/a are in agreement (within errors) with our solution. We used 
our results to reproduce the expected peak frequency Umax and frequency splitting Au, 
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using the formulae of 


Kieldsen fc Bebbind (119951) : 




M/Mg 


(/i/fle)Vr.///8777 


3104 /iHz, 


Ai/ = 

with the reference values from 


M 




1/2 


k 


Mosser et al 


ft 


- 3/2 


138.8 /iHz, 


( 1 ) 


( 2 ) 


(120131) . also used in G13. Our results predict 


Aii = 8.29(20) and 8.69(15) /iHz for the primary and secondary (oscillator in G13 and 
G14), respectively. The peak frequency is however dependent on the temperature, and for 
the discussed range of 4700-5300 K we get 104-98 /iHz for the primary, and 110-104 /iHz 
for the secondary. In particular, the peak frequencies differ by less than the separation Az/ 
and perhaps two modulations have been mixed in the analysis of G13 and G14. 


If we take that only one oscillation is really present, we can hnd the secondary’s 
temperature from the equation [H To get the value of Vmax reported in G14, the secondary 
should have = 5040(260) K, so metallicity close to solar would be preferred 

{[M/H] > —0.39 from Fig. |3]). 


6. Summary 

We have obtained precise mass and radius measurements of two giant components 
of a detached eclipsing binary KIG 09246715. The secondary component was previously 
reported to show solar-like oscillations, and our results conhrm mass and radius derived 
from asteroseismology. We have, however, found that the other star - our primary - has 
very similar properties, thus we suspect that both components may in fact show very 
similar oscillations. We may conclude that the asteroseismology is a promising method of 
measuring stellar masses and radii, and that the correct estimation of these parameters 
for both components of a binary can not be done from light curves only, but still requires 
spectroscopic information. 
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KIC 09246715 



-0.2 -0.1 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 

Phase 


Fig. 1.— The observed (red) and model (blue) Q0-Q17 light curves of K0924. Phase 0 is for 
the deeper eclipse mid-time. Residuals are shown in the lower panel. 
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-0.2 0 0.2 0.4 0.6 

Phase 


Fig. 2.— Radial velocity curves and of KIC 09246715. Filled red circles refer to the primary, 
and open ones to the secondary. The blue lines are the best-htting model curves. Systemic 
velocity is marked by the dotted line. Residuals are shown in the lower panel. Ephemeris 
are the same as in Fig. [TJ 
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Fig. 3.— Left: Comparison of our results with PARSEC isochrones for [M/H] = 0.0 and 
age 950 Myr (red), [M/H] = —0.39 and age of 780 Myr (blue), and [M/H] = 0.20 and 
age 1.04 Gyr (green) on the mass-radius plane. For a given metallicity, the age is restricted 
mainly by the precision in masses. Right: Same isochrones on the temperature-radius plane. 
Grey stripes mark the la ranges of radii we obtained in our analysis. High precision in R 
allows us to estimate temperatures expected for each metallicity/age. With independent 
[M/H] or Teff measurement, coming from spectral analysis for example, the age-metallicity 
degeneracy would be solved. 















Table 1: Measured radial velocities of K0924 with their errors and residuals of the orbital £t 


(all in km s ^). SNR at 5800 A for each observation is also given. 


BJD-2450000 

RVi 

ei (O-C), 

RV 2 

€2 {0-C)2 

SNR 

6865.077724 

38.818 

0.032 

0.022 

-48.417 

0.039 

0.040 

105 

6867.078865 

36.997 

0.036 

0.009 

-46.627 

0.043 

0.001 

75 

6868.138442 

35.671 

0.038 

-0.022 

-45.346 

0.042 

-0.028 

63 

6869.109798 

34.299 

0.043 

-0.023 

-43.941 

0.035 

-0.010 

90 

6914.138901 

-23.176 

0.054 

0.102 

14.398 

0.032 

0.027 

101 

6946.075069 

-25.749 

0.049 

-0.010 

16.872 

0.041 

0.010 

103 

7112.264233 

-26.269 

0.048 

0.005 

17.417 

0.042 

0.013 

81 

7143.230688 

-18.592 

0.038 

-0.049 

9.454 

0.054 

-0.124 

86 
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Table 2: Absolute orbital and physical parameters of K0924. Errors include systematics. 


Parameter 

Value 

± 

P(d) 

171.2770 

0.0006 

To (JD-2454900) 

99.2536 

0.0031 

Tp (JD-2454900) 

81.853 

0.058 

Ki (km s“^) 

33.18 

0.16 

K 2 (km s“^) 

33.58 

0.14 

7i (km s“^) 

-4.643 

0.071 

72 - 7i (km s“^) 

0.15 

0.12 

q 

0.9880 

0.0063 

e 

0.3587 

0.0009 

un 

19.84 

0.44 

ri 

0.04008 

0.00058 

r2 

0.03870 

0.00040 

t (°) 

87.049 

0.031 

J 

1.042 

0.049 

L 2 /L, 

0.964 

0.048 

Ml (Mo) 

2.169 

0.024 

M2 (Mo) 

2.143 

0.025 

Pi (Ro) 

8.47 

0.13 

R 2 (R©) 

8.18 

0.09 

a (Ro) 

211.29 

0.77 

log(^i) 

2.919 

0.013 

log (^2) 

2.944 

0.009 

rmspvi (m s“^) 

45 


rmsRV 2 (m s“^) 

52 


rmspc (mmag) 

0.61 






